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The ^^Fe-specific phonon density of states of Ba(Fei-ccCocc)2As2 single crystals (x = 0.0,0.08) 
was measured at cryogenic temperatures and at high pressures with nuclear-resonant inelastic x- 
ray scattering. Measurements were conducted for two different orientations of the single crystals, 
yielding the orientation-projected ^^Fe-phonon density of states (DOS) for phonon polarizations 
in-plane and out-of-plane with respect to the basal plane of the crystal structure. In the tetragonal 
phase at 300 K, a clear stiffening was observed upon doping with Co. Increasing pressure to 4GPa 
caused a marked increase of phonon frequencies, with the doped material still stiffer than the 
parent compound. Upon cooling, both the doped and undoped samples showed a stiffening, and 
the parent compound exhibited a discontinuity across the magnetic and structural phase transition. 
These findings are generally compatible with the changes in volume of the system upon doping, 
increasing pressure, or increasing temperature, but an extra softening of high-energy modes occurs 
with increasing temperature. First-principles computations of the phonon DOS were performed 
and showed an overall agreement with the experimental results, but underestimate the Griineisen 
parameter. This discrepancy is explained in terms of a magnetic Griineisen parameter, causing an 
extra phonon stiffening as magnetism is suppressed under pressure. 



PACS numbers: 

I. INTRODUCTION 



The recent discovery of superconductivity in iron pnic- 
tides at temperatures up to 55 K has generated great ex- 
citement, and many investigations have sought to clar- 
ify the mechanism responsible for superconductivity in 
these compounds. Although this new class of compounds 
presents some similarities with the cuprate superconduc- 
tors {e.g. layered structure and nearness to magnetism), 
the parent compounds in the case of the iron pnictides 
are metals, rather than insulating oxides. Early first- 
principles calculations of the electron-phonon coupling 
have shown that a conventional electron-phonon coupling 
in the Bardeen-Cooper-Schrieffer theory is too weak to 
explain the magnitude of the observed superconducting 
temperature, Tc [1 . Measurements of phonons in sev- 
eral compounds have been performed, but no significant 
changes in the phonon density of states (DOS) have been 



reported across Tc [2, 3 . However, an isotope effect was 
observed in Bai_ajKa;Fe2As2 and SmFeAsOi-a^Fa^ by Liu 
et a/.[4 , while Shirage et al. observed an inverse isotope 
effect in Bai_a^Ka;Fe2As2 [5j. Also, some anomalies have 
been in reported in the behavior of phonons [6 , and re- 
cent reports have emphasized the potential importance 
of coupling of phonons to the magnetic structure, or to 
spin fluctuations iH H [H HOI [H [H] . 

Superconductivity in the iron pnictides can be achieved 
by doping the parent compounds either in the Fe planes 
themselves, as in the present case of Ba(Fei_a;Coa;)2As2 
[13], or in the separating layers, as in Bai_a;Ka^Fe2As2 or 
LaOi-a^Fa^FeAs. Superconductivity can also be achieved 
in the undoped parent compounds by increasing pres- 
sure [141 llSl lini Hi] 5 although the degree of hydro- 
staticity may influence the appearance of superconduc- 
tivity Several studies have shown that the mag- 
netism in BaFe2As2 is suppressed with increasing pres- 
sure [171 EH [Ej 7 cind that the parent compound becomes 
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superconducting, with a maximum in Tc at P :^ 4GPa 
under quasi-hydrostatic conditions (as obtained in dia- 
mond anvil cehs) [16l [17] . 

Changes in composition, temperature, and pressure 
can all lead to shifts in the phonon energies of a crys- 
tal, but these are often interrelated. For example, tem- 
perature will affect volume through thermal expansion. 
Doping also commonly leads to changes in lattice param- 
eters, besides changing e.g. carrier concentration. As a 
result, the effects of temperature alone (at constant vol- 
ume), or doping alone, is not straightforward to estab- 
lish, unless one records the dependency on the multiple 
coupled parameters. Using Nuclear Resonant Inelastic 
X-ray Scattering (NRIXS), we have performed a system- 
atic study of phonons in Ba(Fei_a;Coa;)2As2, as function 
of temperature, pressure, and doping, in order to com- 
pare the effect of different thermodynamic parameters on 
the phonon DOS. 

Previous investigations have not reported a systematic 
effect of doping on the phonons in Ba(Fei_2;Coa;)2As2, al- 
though effects of doping on phonons have been reported 
in other iron-pnictides [20l ETJ |22]. We show in the 
present study that doping in BaFe2As2 also has an affect 
on the phonons. The effect of an increase in pressure is 
directly observed in the phonons, and the effect of com- 
pression has some similarity with the effect seen upon 
doping. The temperature dependence of the phonons 
can also be interpreted in terms of the change in the 
lattice volume due to thermal expansion, and we could 
compare the predicted effect of thermal expansion within 
the quasi-harmonic approximation, to the effect observed 
upon compression. 



II. SAMPLE SYNTHESIS 

Single-crystalline samples of BaFe2As2 and optimally- 
doped (for maximum Tc) Ba(Fei_a;Co2;)2As2 {x = 0.08) 
were synthesized via a flux-growth technique [13 , using 
an Fe-As flux enriched in the ^^Fe isotope. The isotopic 
enrichment was about 50%. The crystals formed small 
shiny, smooth, plates about 2 mm wide and 0.1 mm thick. 

Resistivity and magnetic susceptibility measurements 
were performed on the ^^Fe-enriched samples. Results 
are shown in Fig. [l] The spin-density wave transition 
temperature of the parent compound was Tsdw = 130 K, 
and the superconducting transition temperature in the 
doped compound was Tc = 22 K. There was no detectable 
shift in the superconducting transition temperature with 
partial ^^Fe enrichment in our samples. Liu et al. have 
reported a positive isotope effect for ^^Fe enrichment in 
Bai_a;Ka;Fe2As2 [4 . However, Shirage et al. reported an 
inverse isotope effect in the same compound with ^^Fe 
and ^^Fe isotopes [5], so a definite isotope effect remains 
to be established. 

X-ray diffraction patterns (A = 1.54056 A) were ac- 
quired at 300 K on powders obtained by crushing the 
platelets. The x-ray patterns were refined with a Rietveld 
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FIG. 1: (Color online) Resistivity curves for BaFe2As2 (122) 
and BaFei.84Coo.i6As2 (Co-122) synthesized with 50% ^'^Fe- 
enriched iron. The top inset shows the derivative of the re- 
sistivity, and the magnetization, for the parent compound 
BaFe2As2. The bottom inset shows the resistivity and mag- 
netization of BaFei.84Coo.i6As2 at low temperature. 



TABLE I: Rietveld refinements of crystal structure of Fe- 
enriched samples (T = 300 K, A = 1.54056A). 



sample 


BaFe2As2 


BaFei.84Coo.i6As2 


a = 6 (A) 


3.9613(1) 


3.9604(1) 


c(A) 


13.0148(3) 


12.9793(3) 


* 


0.35423(7) 


0.35423(8) 



* Fractional coordinate. 



procedure, as implemented in the software FULLPROF 
[23 . Results are listed in Table |lj The substitution of 
Fe by Co leads to a slight contraction of c parameter 
(—0.27%), while the a and b parameters are almost un- 
changed (difference smaller than 0.03%). The contrac- 
tion of the c-axis upon doping and the lack of change for 
a and b axes are consistent with the lattice parameters 
previously reported for the case of samples synthesized 
from natural iron [13]. 



III. NUCLEAR-RESONANT INELASTIC X-RAY 
SCATTERING MEASUREMENTS 

Nuclear resonant inelastic x-ray scattering (NRIXS) 
measurements [lH |25l [26] were performed at cryogenic 
temperatures at beamline 3-IDD and high pressures at 
beamline 16-IDD (HP-CAT) at the Advanced Photon 
Source at the Argonne National Laboratory. 

In all measurements, the incident photon energy was 
tuned to 14.4124 keV, the nuclear resonance energy of 
^^Fe. The NRIXS signal was measured with multiple 
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FIG. 2: (Color online) NRIXS data measured on BaFe2As2 
sample at T = 300 K and P = OGPa (solid line, "sample"), 
and measured resolution function of beamline 16-IDD (dashed 
line, "resolution"). AE is the shift from the ^^Fe nuclear 
resonance energy. 



avalanche photodiode detectors positioned 90° from the 
direction of the beam, with the exception of the IP high 
pressure measurement where detectors were positioned 
60° from the direction of the beam (due to pressure cell 
geometry constraints). Data were collected in scans of in- 
cident photon energy, with = —60 to +60 meV from 
the resonant energy, in steps of 0.5 meV. The experimen- 
tal energy resolution function was measured with a single 
avalanche photodiode placed in the forward beam direc- 
tion, recording the intensity as a function of the shift of 
the incident energy away from the ^^Fe resonance energy 
(the data for the instrument resolution were summed over 
all runs performed in the same conditions). Representa- 
tive spectra for the NRIXS data and the instrumental 
resolution are shown in Fig. |2] The monochromator en- 
ergy resolution (full width at half maximum) was 2.2 meV 
in the measurements as function of pressure, and 2.5 meV 
in the measurements as function of temperature. 

All of the NRIXS data reduction were performed us- 
ing the standard software PHOENIX [27,. The raw 
NRIXS spectra, given as intensity versus the angle of 
the monochromator crystals, were converted to intensity 
versus energy transfer. The first few bins on the low en- 
ergy side (-60meV to -55meV) were used to determine 
an energy independent background, and was removed for 
all energy transfers. The elastic peak was removed using 
the measured resolution function. The contribution from 
multiphonon scattering processes was subtracted using a 
self-consistent procedure based on a Fourier-log method 
|28] , and the Fe-partial phonon density of states was ob- 
tained by correcting for the thermal occupation factor. 



A. Pressure-Dependent Measurements 

Measurements as function of pressure were performed 
in panoramic diamond-anvil cells (DAG), with diamonds 



of 500 /im-diameter culets. The samples were 100 /urn- 
wide plate-like crystals, and were contained inside the 
DAC's by Be gaskets drilled with 130 jam holes (platelets 
remained parallel to the gasket plane). The pressure 
medium was silicone oil, which provides a nearly hydro- 
static pressure medium up to 10 GPa [29]. The pressure 
inside the DAC's was determined through the fluores- 
cence line of ruby crystals loaded with the samples in 
the pressure medium [30] . Measurements were performed 
with the cells at ambient pressure (0 GPa) and with the 
cells pressurized to 4.0 ± 0.2 GPa. 
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FIG. 3: Schematic of geometry of high pressure NRIXS mea- 
surements. See text for description. 

Out-of-plane (OP) measurements were performed with 
the incident beam passing through the diamonds and the 
scattered NRIXS signal measured through the Be gasket. 
In-plane (IP) measurements were performed with the in- 
cident beam passing through the Be gasket and the scat- 
tered NRIXS signal also measured through the Be gasket. 
The geometry of the high-pressure NRIXS measurements 
is illustrated in Fig. [3j 



B. Temperature-Dependent Measurements 

Measurements as function of temperature (10 K to 
300 K, at ambient pressure) were performed in an evac- 
uated He-flow cryostat, with the sample mounted on a 
Gu cold finger and exposed to the x-ray beam through a 
beryllium dome. The temperature was monitored with 
two thermocouples positioned respectively on the cryo- 
stat head, and directly next to the sample on a Gu 
bracket. The readings from the two thermocouples were 
in good agreement throughout the measurements. 
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IV. EXPERIMENTAL RESULTS 
A. Effect of Doping 

The Fe-partial phonon DOS (300 K, OGPa) for IP and 
OP modes in BaFe2 As2 and BaFei.84Coo.i6As2 are shown 
in Fig. |4] The Fe-partial phonon DOS is composed of 
three peaks around 13, 25, and 32meV. The low and 
medium-energy peaks are more prominent in the OP- 
polarized DOS, indicating that a larger proportion of 
displacements along the c axis in these energy ranges. 
This is physically intuitive, as the crystal structure is ex- 
pected to have lower vibration frequencies for displace- 
ments of Fe atoms along c. In the parent compound 
(x = 0) , we find for the first-moment of the phonon ener- 
gies: (Eff^) = 24.57 ± 0.04 meV for IP vibrations, while 
i^of) = 22.93 ± 0.06 meV for OP vibrations. 

As can be seen in this Fig. |4j the phonons are stiffer 
in the doped material than in the parent compound, for 
both the IP and OP-polarized vibrations. This is par- 
ticularly well seen on the peak at about 32 meV and at 
the phonon cutoff. The mass difference of Co and Fe 
atoms cannot account for this shift, since Co is more 
heavy than Fe. The phonon stiffening is about 0.5 meV 
at the cutoff ( ~ 1.5% change), which is significant, con- 
sidering that doping only causes a modest contraction 
in the c lattice parameter of order 0.25%, and almost 
no change in the a and b parameters. The stiffening 
of phonons upon doping was also observed in the mea- 
surements at 4GPa, where a similar magnitude in the 
shift upon doping persisted. Our analysis of the cut- 
off energy gives AEl^^ = 0.5 ± 0.1 meV for IP modes, 
and A£;^P = 0.4 ± 0.1 meV for OP modes, at 300 K and 
OGPa. 



B. Effect of Pressure 

Measurements under pressure were performed with sin- 
gle crystals loaded inside the diamond anvil cells. The 
single crystals were oriented with the c-axis parallel to 
the axis of the DAC. This orientation was verified using 
in-situ diffraction (both at CPA and 4 GPa) , and also 
confirmed by comparing the phonon spectra at GPa 
with samples in the DACs and the measurements with- 
out the pressure cells. Area diffraction patterns indi- 
cated some turbostratic disorder in the (a, b) plane of 
the samples (rotational stacking disorder along c), but 
this type of disorder does not hinder the measurement 
of separate in-plane and out-of-plane polarizations with 
our setup. The results for the Fe-partial phonon DOS 
of BaFei.84Coo.i6As2 and BaFe2As2 are shown in Figs.jS] 
and[6j respectively. 

The effect of pressure on the phonon DOS is signifi- 
cant. In both the parent and the doped samples, a large 
stiffening is observed. The pressure dependence of the in- 
plane (IP) modes is similar in doped and undoped ma- 
terials, but the out-of-plane modes (OP) appear more 
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FIG. 4: (Color online) (a) Fe-partial phonon DOS for modes 
with in-plane (IP) polarizations in Ba(Fei-xCocc)2As2 for x = 
0.0 and x = 0.08, at ambient conditions (T = 300 K, P = 
OGPa). (b) Same for out-of-plane (OP) polarizations. 
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FIG. 5: (a) Partial Fe phonon DOS for in-plane (IP) polar- 
ized modes in BaFei.84Coo.i6As2, at GPa and 4.0 di 0.2 GPa 
(300 K). (b) Same for out-of-plane polarizations (OP). 

sensitive to pressure in the doped material. This may be 
due to the doping-induced c-axis reduction in the super- 
conductor. From the slopes, we find the pressure coef- 
ficients: d{Eop)/dP = 0.68 and 0.45meV/GPa for OP 
modes in doped and parent compounds, respectively, and 
d{Eip)/dP = 0.40 and 0.43meV/GPa for IP modes in 
doped and parent compounds, respectively. 

Our in-situ diffraction data on BaFei.84Coo.i6As2 and 
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FIG. 6: (a) Partial Fe phonon DOS for in-plane (IP) polarized 
modes in BaFeaAsa, at OGPa and 4.0 =b 0.2 GPa (300 K). (b) 
Same for out-of-plane polarizations (OP). 
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FIG. 7: (Color online) (a) Average energy as function of pres- 
sure in BaFe2As2 and BaFei.84Coo.i6As2 (300 K). (b) Energy 
of top optical peak as function of pressure in BaFe2As2 and 
BaFei.84Coo.i6As2 (300 K). 



BaFe2 As2 at 4 GPa show lattice compressions compared 
to GPa that are in good agreement with those reported 
by Kimber et al. [31 (although we could not resolve 
in-situ the difference in c-parameter between the two 
compounds). We find the following lattice parameters 
for BaFei.84Coo.i6As2 at 4 GPa and 300 K: a = b = 
3.90 ± 0.01 A ; c = 12.71 ± 0.01 A. Within the error 
bars of our in-situ measurements, the same values are 



TABLE II: Experimental Griineisen parameters of ^^Fe- 
enriched Ba(Fei_xCocc)2As2, obtained from NRIXS measure- 



ments as function 


of pressure (T= 


300 K) or 


temperature 


(P= OGPa). 










0'T=300 k 


^X=300 K 


^P = GPa 


polarization 


IP 


OP 


OP 


{E) 








BaFe2As2 


1.4 ±0.2 


1.6 ±0.2 


1.5 ±0.3 


BaFe1.8Coo.2As2 


1.3 ±0.2 


2.1 ±0.2 


1.9 ±0.6 


^optic,2 (32meV) 








BaFe2As2 


1.8 ±0.1 


1.7±0.1 


2.2 ±0.2 


BaFe1.8Coo.2As2 


1.7±0.1 


1.8 ±0.1 


2.9 ±0.3 



obtained for BaFe2As2 at 4GPa and 300K. Both the su- 
perconductor and the parent compound appear to com- 
press by the same relative amount between OGPa and 
4 GPa: -2.5 ± 0.1% for c, and about -1.5 ± 0.2% for a, b 
(all at 300K). Kimber et al. have reported that the com- 
pression of a and c axes is very linear with pressure from 
1 GPa to 6 GPa, and that compressibilities along both 
axes are unchanged between 17 and 150 K. Our results 
indicate that this also holds at 300 K. 

From the experimental change in volume, and the 
measured phonon frequencies, we obtained the mode- 
dependent Griineisen parameters at T = 300 K, 7 = 
—d\nE/d\nV^ for the average phonon energy as well as 
the high-energy optical modes. Results are listed in Ta- 
ble [HI 



C. Effect of Temperature 

The OP Fe-partial phonon DOS of BaFei.84Coo.i6As2 
was measured at 15, 40, 160, and 300 K, and that of 
BaFe2As2 was measured at 10, 100, 160, and 300 K. We 
do not observe any large change in the OP Fe-pDOS of 
BaFei.84Coo.i6As2 across Tc = 22 K. The main effect ob- 
served is an overall softening of the phonon DOS with 
increasing T. This softening corresponds to a downward 
shift of about 1 meV in the position of the 32 meV peak 
between 10 K and 300 K, compatible with the effect ob- 
served on the same peak in LaOo.gFo.iFeAs [2 . As we 
show below, this behavior is compatible with the quasi- 
harmonic (QH) model, according to which phonon ener- 
gies depend on temperature through the thermal expan- 
sion: AE{T)/Eo = -7 AV(T)/Fo, with 7 the Griineisen 
parameter. 

The parent compound also shows a systematic de- 
crease in energy of all the modes upon heating from 
10 K to 300 K, by an amount comparable to that in 
the doped material. The phonon DOS in the tetrago- 
nal and orthorombic phases are similar overall, as ex- 
pected from the small magnitude of the structural dis- 
tortion across the transition [32] |33]. However, we do 
observe a slight change in shape of the DOS across the 
tetragonal-orthorombic phase transition at 130 K. The 
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FIG. 8: (Color online) (a) OP Fe-pDOS of BaFei.84Coo.i6As2 
as function of temperature, as measured with NRIXS. (b) 
Same for BaFe2As2. The curves for measurements at 160 and 
300 K in the lower panel were vertically offset for clarity. 



optical phonon peak at ~ 32 meV appears to be the most 
sensitive, and shifts down abruptly across the transition 
(heating). Additionally, the phonon peak at 12meV ap- 
pears more broad above the phase transition. 

The temperature dependence of the average energy, 
and of the high-energy optical phonon (around 32 meV) 
for OP modes is plotted in Fig. |9] The average energy 
{E) was obtained from the first moment of the measured 
DOS, and the high-energy optical peak position was es- 
timated by fitting the top half of the peak to a Gaus- 
sian. Error bars are from experimental counting statis- 
tics. As can be seen on this figure, the phonon energy 
increases with decreasing temperature in both BaFe2As2 
and BaFei.84Coo.i6As2, as expected from the decrease in 
overall volume of the lattice at lower temperatures. The 
behavior in BaFei.84Coo.i6As2 is smooth, and follows the 
behavior expected from QH theory. Our data show the 
possibility of a small softening upon cooling across Tc (see 
Fig. [9|, compatible with the upturn in volume expansion 
reported in [34", 'SF. In the case of BaFe2As2, there is a 
discontinuity in the frequency with a jump at the tem- 
perature of the tetragonal-orthorombic (and magnetic) 
transition. This jump is 0.8 ±0.2 meV for the high-energy 
optic OP mode, with the orthorombic phase stiffer than 
the tetragonal phase. Our observation can be related to 
the smaller volume of the orthorombic phase: Bud'Ko et 
al have reported a volume collapse across Tsdw , mostly 
associated to a decrease in the basal plane dimensions 
[34] . When both the doped and undoped materials are 
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FIG. 9: (Color online) (a) Temperature dependence of average 
energy for Fe vibrations with out-of-plane polarization (OP), 
in BaFe2As2 and BaFei.84Coo.i6As2. (b) Same for OP high- 
energy optic peak. Markers with error bars are experimental 
values. Dashed curves give the calculated T-dependence for 
a quasiharmonic (QH) model, ^qh(T) = Eo{V{T)/Voy^ , 
using thermal expansion data reported in [MiES], see text. 



in the tetragonal phase, the high-energy optic mode in 
the doped material is stiffer than in the undoped com- 
pound by about 0.5 meV. Below Tsdw, the orthorombic 
BaFe2As2 is stiffer than tetragonal BaFei.84Coo.i6As2. 

We note that Choi et al. have reported Raman mea- 
surements in CaFe2As2, which show a 0.5 meV step-like 
stiffening of the Big Raman mode (at 26 meV) when cool- 
ing across the structural transition in that compound 
j36] . Our data around the energy of the Big Raman mode 
(26meV) in BaFe2As2 are less clear than for the 31meV 
mode, but show an energy difference of 0.3 ±0.2 meV be- 
tween 160 K and 100 K. The measurements of Chauviere 
et al. for the Big Raman mode in BaFe2As2 do not 
show a clear discontinuity either, but seem to show a 
kink around Tsdw [37 . We do not observe any trace 
of the splitting of the low-energy Eg Raman mode re- 
ported by these authors [37 , but this mode corresponds 
to vibrations in the (a, 6) plane (IP-polarized) [3F. 

The temperature dependence of the phonon frequen- 
cies can be computed from the thermal expansion of the 
crystal using the quasi-harmonic approximation. Using 
the temperature dependence of the cell volume reported 
by Bud'Ko et al. for the doped compound [34], and by 
Daluz et al. for the parent compound [35 , we determined 
the Griineisen parameters giving the best fits to our mea- 
surements of the 32 meV OP optical mode (corresponding 
to dashed curves in Fig. [9| . Assuming a T-independent 
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Griineisen parameter 7 = —d\nE/d\nV = const., we 
optimized 7 and the energy at low temperature, Eo{T = 
10 K), so that Equ{T) = Eo{V{T)/Voy^ would min- 
imize the of ^qh(^) — ^NRixs(^) for all measured 
temperatures. Results are listed in Table |TTj We find 
7|p=o = 2.2 ±0.2 for BaFe2As2 and 7|p=o = 2.9 ±0.3 for 
BaFei.84Coo.i6As2 for the high-energy optical peak with 
OP polarization. It is to be noted that these "thermal" 
Griineisen parameters are significantly larger than the 
pure volume Griineisen parameters obtained from mea- 
surements under pressure (7|t=300 — 1-7 for the high- 
energy OP peak), or from the values obtained from DFT 
calculations, as discussed below. 



V. FIRST-PRINCIPLES SIMULATIONS 

The first principle calculations of in-plane (IP) and 
out-of-plane (OP) Fe-partial phonon DOS for BaFe2As2 
were done within the framework of density functional 
theory (DFT), with the generalized gradient approxima- 
tion of Perdew et al. [39 The phonon calculations were 
done using linear response as implemented in Quantum 
Espresso [40 . We used experimental lattice parameters 
when available but relaxed the As height via energy min- 
imization. The basis set cutoff for the wave functions 
was 40 Ry, while a 400 Ry cutoff was used for the charge 
density. A 4 x 4 x 4 special /c-point grid was used for 
Brillouin zone sampling. The dynamical matrices were 
calculated on a 4 x 4 x 4 grid and the dynamics was then 
interpolated on a much finer 48 x 48 x 48 grid to obtain 
the phonon DOS. 

The lattice parameters used in the calculation were: 
a = b = 3.9625 A, c = 13.0168 A, corresponding to 
the experimental structure at ambient conditions (300 K, 
GPa) [32 . Relaxation of the As fractional coordi- 
nate along c gave zas = 0.34515. For the comparison 
with measurements at 4 GPa, computations were per- 
formed on a compressed unit cell with a = b = 3.913 A, 
c = 12.698 A, corresponding the high-pressure diffrac- 
tion measurements reported in [31 , and also in good 
agreement with our experimental estimates. In this case, 
the relaxed zas was 0.34786. In order to compare the 
phonon DOS of the parent compound BaFe2As2 and the 
doped compound BaFe1.8Coo.2As2, additional computa- 
tions were performed on a unit cell of BaFe2As2 com- 
pressed along the c axis by 0.25%. The slight compres- 
sion of the c-axis induced almost no change in the relaxed 
fractional As position, zas = 0.34514. The purpose of 
this calculation was to study the effect of geometrical 
changes of the structure on the phonons, without chang- 
ing the composition. 

Results for the Fe-partial DOS for in-plane and out- 
of-plane polarizations at GPa and 4 GPa are shown in 
Fig. [lO| The computed phonon DOS curves were con- 
volved with the measured experimental resolution func- 
tion to compare with the experimental DOS. The com- 
puted phonon DOS curves are in good agreement with 
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FIG. 10: (a) Fe-partial DOS computed with DFT, for in-plane 
(IP) polarizations in BaFe2 As2 at GPa and 4 GPa, convolved 
with the experimental resolution function, (b) Same for out- 
of-plane polarizations (OP). 



the measurements, for both the shape of the DOS and 
the energies of the modes. The number of peaks and their 
relative intensities agree well with our experimental data. 
The DFT frequencies are consistently about 10% stiffer 
than the measured frequencies, however. This short- 
coming is well known for non-spin-polarized calculations 
[3 [3 [H]. A treatment with spin-polarized wavefunc- 
tions typically brings better agreement with measured 
phonon frequencies [7, 8 , however the agreement within 
the non-spin-polarized calculations is sufficient to inves- 
tigate trends as function of unit cell volume and geom- 
etry. From the average energy of the computed phonon 

(E) 

DOS, we find the Griineisen parameters 7jp = 1.11 and 
7ot ~ ^'^0, while for the high-energy optic modes, we 
obtain 7^p*^^ = 1.58 and 7op^^^ — l-^^- These theoretical 
values are systematically lower than our experimental re- 
sults (see Table [TT|) . This discrepancy indicates that the 
suppression in magnetism upon compression may result 
in an additional stiffening of phonons, which is not cap- 
tured in the present computations. 



The projected phonon DOS computed on BaFe2As2 
compressed along the c axis is compared to the computed 
DOS at the experimental (uncompressed) geometry in 
Fig. [11] As can be seen in this figure, the compression 
of the c-axis alone leads to a stiffening of phonon modes 
polarized both along c (OP) and perpendicular to it (IP), 
as observed experimentally upon Co-doping. The mag- 
nitude of the stiffening for both polarizations is also in 
good agreement with the experimental observation. 
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FIG. 11: (a) Fe-partial DOS computed with DFT, for in-plane 
(IP) polarizations in BaFe2As2 at 300 K with experimental 
unit cell, and for BaFe2As2 with cell compressed 0.25 % along 
c axis. Curves are convolved with the experimental resolution 
function, (b) Same for out-of-plane (OP) polarizations. 



VI. DISCUSSION 

Our phonon measurements as function of doping, pres- 
sure and temperature show a systematic trend. In 
all cases, we observe a strong correlation between the 
phonon energies and the unit cell volume, with contrac- 
tions in the lattice parameters leading to stiffer phonons. 
This behavior is expected on the basis of inter-atomic 
force-constants weakening with increasing atomic sepa- 
ration. The change in phonon frequencies is related to 
the change in volume through the Griineisen parameter 
(quasiharmonic approximation), and we have determined 
the Griineisen parameters of IP and OP polarizations, for 
both the average phonon energy and the higher-energy 
optical peak. The Griineisen parameters obtained from 
measurements under pressure (listed in Table [ll| repre- 
sent the dependence of phonon energy levels on pressure 
alone (at 300 K). On the other hand, the Griineisen pa- 
rameters derived from our T-dependent measurements 
and the reported thermal expansion combine the effect 
of the volume change and other T-induced phenomena, 
such as carrier excitation across the Fermi level, or T- 
dependent magnetic excitations. 

Upon closer examination, the contraction of the unit 
cell with Co doping (0.25% in volume) at 300 K is ex- 
pected to result in a stiffening of the average phonon 
energy of about 0.4%, A{E)/{E) - 0.004, based on 
our measured Griineisen parameters. The observed 
stiffening is in good agreement with this result, with 
A{Eop)/{Eop) = 0.004 ±0.002 (beamline 16-IDD) and 
0.006 ± 0.002 (beamline 3-IDD), and A(£^ip)/(£^ip) = 
0.01 ±0.004 (beamline 16-IDD). Thus, we conclude that 
electron doping results in a stiffening of phonons that 



can be accounted for based on the volume contraction 
alone. This is in agreement with the results of the DFT 
simulations, which show that a 0.25% compression in the 
c lattice parameter accounts for the observed change of 
both the OP and IP-polarized phonon DOS upon doping 
(see Fig. 11 ). In both theory and experiment, the peak at 
32 meV appears to be the most sensitive to compressions 
of the c axis (for IP as well as OP polarizations). 

Recently, Kimber et al [31] have reported a similarity 
between the effects of pressure and doping on the lattice 
constants of (Ba,K)Fe2As2. Our measurements indicate 
that doping and compression also have similar effects on 
the phonons, which is easily understood as Co exerting 
an effective chemical pressure (Co is a smaller species 
than Fe), which leads to a smaller unit cell and stiffer 
phonons. 

The Griineisen parameters obtained from first- 
principles calculations (corresponding to lattice compres- 
sions) are significantly smaller than the experimental val- 
ues determined from high-pressure phonon data. Al- 
though the non-spin-polarized calculations are expected 
to predict stiffer phonon frequencies than the experimen- 
tal values, this should be the case at all volumes, and does 
not account for the underestimated Griineisen parame- 
ter. A possible explanation is that compression leads to 
a weakening of the magnetic moments in the material, 
and an extra phonon stiffening, compared to a pure vol- 
ume effect. Since our calculations are non- magnetic at all 
volumes, this extra stiffening would be missed, leading to 
smaller predictions for the Griineisen parameters. 

Using a recently formulated Landau theory [10], this 
effect can be estimated for the case of the As Raman 
phonon mode, which corresponds to As planes vibrating 
against the Fe planes. We may write the magnetic free 
energy as: 



Fm = AM'^ ± BM^ - {a{z - Zc) - P{z - Zcf)M'^ 

^^{z-{z)f , (1) 

where M is the Fe moment, z is the separation between 
Fe and As planes, Zc is the Stoner quantum critical point 
in z for onset of magnetism, (z) is the thermal average 
of and Kq is the (harmonic) elastic constant without 
the spin-lattice interaction. By minimizing Eq. 1 (^4=0 
at 2:c) , we obtain: 



M 



I a 



{z- Zc) - (3{z- Zcf 



2B 



1/2 



(2) 



which agrees with calculation [10 and experiment [42]. 
Also, by expanding Eq. 1 in z we obtain the elastic con- 
stant renormalized for the spin-lattice interaction. 



(3) 



Using the values evaluated by the DFT calculation 
for Ba(Feo.92Coo.o8)2As2 = 1.278 A, aj^B = 
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11.67 /i|/A, a = 0.193 eV/A/i|, (3 = 0.137eV/A /i| , 
and {z) = 1.36 A, we obtain the enhancement to the 
Griineisen constant due to the loss of magnetism under 
increasing pressure, 



A7 = 



IdlnK 
2 dhiz 



1.81. 



(4) 



The difference in the Griineisen parameter obtained from 
DFT calculations without the spin-lattice effect and the 
values measured by NRIXS for BaFei.84Coo.i6As2 is 
A7 ~ 0.8, about a factor of two smaller than the value de- 
rived above. However, this overestimate is expected since 
the measured value of 7 is an average over many phonon 
modes, while the model calculation was for the phonon 
mode which is most strongly affected by spin-phonon cou- 
pling. Also, the Fe-As bond is hard, and compresses less 
than the lattice parameters under pressure [31 . 

The temperature dependence of phonon energies is 
compatible with the increase in volume upon thermal 
expansion, as was shown in Fig. [9) However, the iso- 
baric Griineisen parameter for high-energy optic modes 
(obtained from our data and reported thermal expansion 
curves), 7|p 2.6, (P = GPa) is significantly larger 
than its isothermal (T = 300 K) counterpart, 7|t — 1.6. 
It thus appears that T-induced modifications of the elec- 
tronic and/or magnetic structure may contribute directly 
to the T dependence of phonons. It is possible that 
a weakening of magnetism with increasing temperature 
causes an enhancement in 7, similarly to the case of high 
pressure discussed above. 

Also, since the electronic DOS of BaFe2 As2 has a large 
peak in the DOS below the Fermi level, it is possible that 
increased temperature leads to extra free carriers, yield- 
ing larger screening and softer phonons than would be ob- 
tained through a pure volume expansion. We point out 
that the T-dependence of phonons can exhibit system- 
atic departures from the quasiharmonic model, when the 
Fermi level is in proximity of sharp features in the elec- 
tronic DOS [43, 44 . In particular, metals whose Fermi 
level lies in a pseudo-gap in proximity of a peak in the 
electronic DOS, as is the case in the FeAs compounds [T3^ 
(and also in Cr and Mo) can exhibit an excess phonon 
softening with increasing temperature. The dependence 
of phonons on temperature in Ba(Fei_a;Coaj)2As2 thus 



presents some similarity with that of Cr and Mo, al- 
though this dependence is not altered over the small 
Co doping range which induces the superconductivity in 
Ba(Fei_^Co^)2As2. 



VII. CONCLUSIONS 

Using NRIXS, we have measured the partial phonon 
DOS for in-plane and out-of-plane vibrations of ^''Fe 
atoms in BaFe2As2 and BaFe1.8Coo.2As2, as function of 
temperature and pressure. These experimental phonon 
DOS curves were compared to density functional calcu- 
lations. Good agreement was observed in the trend as 
function of doping and pressure between the experimen- 
tal and theoretical results. The dependence of phonons 
on doping, pressure, and temperature can be understood 
primarily based on the dependence of the lattice con- 
stants on these thermodynamic parameters. However, 
magnetism appears to induce an additional phonon stiff- 
ening upon increasing pressure. Also, the dependence of 
phonons on temperature is more pronounced than would 
be expected for a pure volume thermal expansion, for 
both the parent and doped compounds. This extra soft- 
ening may be attributed to thermally-induced rearrange- 
ments of electronic or magnetic structures. 



VIII. ACKNOWLEDGEMENTS 

We thank Steve Nagler for help with providing the 
isotopically-enriched iron. Experimental work at ORNL 
was supported by the Scientific User Facilities Division 
and the Division of Materials Sciences and Engineer- 
ing, Office of Basic Energy Sciences, DOE. Theory work 
at ORNL was supported by DOE, Division of Mate- 
rials Sciences and Engineering. Part of the research 
was performed by a Clifford G. Shull Fellow (O.D.) and 
by Eugene P. Wigner Fellows (A.S.S. and M.A.M.) at 
ORNL. Use of the HPCAT facility was supported by 
DOE-BES, DOE-NNSA (CDAC), NSF, DOD TACOM, 
and the W.M. Keck Foundation. Use of the APS was 
supported by DOE-BES, under Contract No. DE-AC02- 
06CH11357. 



[1] L. Boeri, O. V. Dolgov and A. A. Golubov, Phys. Rev. 

Lett. 101, 026403 (2008). 
[2] A. D. Christiansen, M. D. Lumsden, O. Delaire, M. B. 

Stone, D. L. Abernathy, M. A. McGuire, A. S. Sefat, R. 

Jin, B. C. Sales, D. Mandrus, E. D. Mun, P. C. Canfield, 

J. Y. Y. Lin, M. Lucas, M. Kresch, J. B. Keith, B. Fultz, 

E. A. Goremychkin, and R. J. McQueeney, Phys. Rev. 

Lett. 101, 157004 (2008). 
[3] D. Phelan, J. N. Millican, E. L. Thomas, J. B. Leao, Y. 

Qiu, and R. Paul, Phys. rev. B 79, 014519 (2009). 
[4] R. H. Liu, T. Wu, G. Wu, H. Chen, X. F. Wang, Y. L. 



Xie, J. J. Ying, Y. J. Yan, Q. J. Li, B. C. Shi, W. S. Chu, 

Z. Y. Wu, and X. H. Chen, Nature 459, 64 (2009). 
[5] P.M. Shirage, K. Kihou, K. Miyazawa, C.-H. Lee, H. 

Kito, H. Eisaki, T. Yanagisawa, Y. Tanaka, A. lyo, 

arXiv:0903.3515Y2. 
[6] R. Mittal, L. Pintschovius, D. Lamago, R. Heid, K-P. 

Bohnen, D. Reznik, S. L. Chaplot, Y. Su, N. Kumar, S. 

K. Dhar, A. Thamizhavel, and Th. Brueckel, Phys. Rev. 

Letters 102, 217001 (2009). 
[7] T. Yildirim, Physica C 469, 425-441 (2009). 
[8] M. Zbiri, H. Schober, M. R. Johnson, S. Rols, R. Mittal, 



10 



Y. Su, M. Rotter, D. Johrendt, Phys. Rev B 79, 064511 
(2009). 

[9] D. Reznik, K. Lokshin, D. C. Mitchell, D. Parshall, W. 
Dmowski, D. Lamago, R. Held, K.-P. Bohnen, A.S. Sefat, 
M. A. McGuire, B. C. Sales, D. G. Mandms, A. Subedi, 
D. J. Singh, A. Alatas, M. H. Upton, A. H. Said, A. 
Cunsolo, Yu. Shvydko, T. Egami, arXiv:0908.4359Y2. 

[10] T. Egami, B. V. Fine, D. J. Singh, D. Parshall, C. de la 
Cruz, P. Dai, arXiv:0908.436lYl. 

[11] S. E. Hahn, Y. Lee, N. Ni, P. C. Canfield, A. I. Gold- 
man, R. J. McQueeney, B. N. Harmon, A. Alatas, B. M. 
Leu, E. E. Alp, D. Y. Chung, L S. Todorov, and M. G. 
Kanatzidis, Phys. Rev. B 79, 220511(R) (2009), 

[12] R. M. Fernandes, L. H. VanBebber, S. Bhattacharya, P. 
Chandra, V. Keppens, D. Mandrus, M. A. McGuire, B. 
C. Sales, A. S. Sefat, J. Schmalian, a rXiv:091 1.3084vl. 

[13] A. S. Sefat, R. Jin, M. A. McGuire, B. C. Sales, D. J. 
Singh, and D. Mandrus, Phys. Rev. Lett. 101, 117004 

(2008) . 

[14] M. S. Torikachvih, S. L. Bud'ko, N. Ni, and P. C. Can- 
field, Phys. Rev. Lett. 101, 057006 (2008) 
[15] T. Park, E. Park, H. Lee, T. Khmczuk, E. D. Bauer, 

F. Ronning, and J. D. Thompson, J. Phys.: Condens. 
Matter 20, 322204 (2008). 

[16] P. L. Alireza, J. Gillett, Y. T. Chris Ko, E. Suchitra, and 

G. G. Lonzarich, J. Phys.: Condens. Matter 21, 012208 

(2009) , 

[17] E. Colombier, S. L. Bud'ko, N. Ni, and P. C. Canfield, 

Phys. Rev. B 79, 224518 (2009). 
[18] K. Matsubayashi, N. Katayama, K. Ohgushi, A. Ya- 

mada, K. Munakata, T. Matsumoto, Yoshiya Uwatoko, 

arXiv:0905.0968v2. 
[19] Hideto Fukazawa, Nao Takeshita, Takehiro Yamazaki, 

Kenji Kondo, Kenji Hirayama, Yoh Kohori, Kiichi 

Miyazawa, Hijiri Kito, Hiroshi Eisaki, and Akira lyo, J. 

Phys. Soc. Jpn. 77, 105004 (2008). 
[20] M. Le Tacon, M. Krisch, A. Bosak, J.-W. G. Bos and S. 

Margadonna, Phys. Rev. B 78, 140505 (R) (2008). 
[21] M. Le Tacon, T. R. Forrest, Ch. Riiegg, A. Bosak, A. C. 

Walters, R. Mittal, H. M. R0nnow, N. D. Zhigadlo, S. 

Katrych, J. Karpinski, J. P. Hill, M. Krisch, D.F. Mc- 

Morrow, arXiv:0909.0913vl. 
[22] Chul-Ho Lee, Kunihiro Kihou, Kazumasa Horigane, 

Satoshi Tsutsui, Tatsuo Fukuda, Hiroshi Eisaki, Akira 

lyo, Hirotaka Yamaguchi, Alfred Q. R. Baron, Markus 

Braden, Kazuyoshi Yamada, |arXiv:0911.3215 vl. 
[23] J. Rodriguez-Carvajal, Physica B 192, 55 (1993). 

http: / / www.ill.eu / sites / fullprof / 
[24] E. E. Alp, W. Sturhahn, T. S. Toellner, J. Zhao, M. 

Hu and D. E. Brown, Hyperfine Interactions 144/145 

(2002). 



[25] M. Seto, Y. Yoda, S. Kikuta, X.W. Zhang, and M. Ando, 

Phys. Rev. Letters 74, 3828 (1995). 
[26] W. Sturhahn, T.S. Toellner, E.E. Alp, X. Zhang, M. 

Ando, Y. Yoda, S. Kikuta, M. Seto, C.W. Kimbah, and 

B. Dabrowski, Phys. Rev. Letters 74, 3832 (1995). 
[27] W. Sturhahn, Hyperfine Interact. 125, 149 (2000). 
[28] D.W. Johnson and J.C.H. Spence, J. Phys. D 7, 771 

(1974). 

[29] S. Klotz, J.-C. Chervin, P. Munsch, and G. Le Marchand, 

J. Phys. D: Appl. Phys. 42 (2009). 
[30] G.J. Piermarini, S. Block, J.D. Barnett, R.A. Forman, 

Journal of Applied Physics 46(6), 2774 (1975) 
[31] S. A. J. Kimber, A. Kreyssig, Y.-Z. Zhang, H. O. Jeschke, 

R. Valenti, F. Yokaichiya, E. Colombier, J. Yan, T. C. 

Hansen, T. Chatterji, R. J. McQueeney, P. C. Canfield, 

A. I. Goldman and D. N. Argyriou, Nature Materials 8, 

471 (2009). 

[32] M. Rotter, M. Tegel, D. Johrendt, I. Schellenberg, W. 
Hermes, and R. Pottgen, Phys. Rev. B 78, 020503(R) 
(2008). 

[33] Q. Huang, Y. Qiu, Wei Bao, M. A. Green, J. W. Lynn, 

Y. C. Gasparovic, T. Wu, G. Wu, and X. H. Chen, Phys. 

Rev. Letters 101, 257003 (2008). 
[34] S. L. Bud'ko, N. Ni, S. Nandi, G. M. Schmiedeshoff, and 

P. C. Canfield, Phys. Rev. B 79, 054525 (2009). 
[35] M. S. da Luz, J. J. Neumeier, R. K. Bollinger, A. S. Sefat, 

M. A. McGuire, R. Jin, B. C. Sales, and D. Mandrus, 

Phys. Rev. B 79, 214505 (2009). 
[36] K.-Y. Choi, D. Wulferding, P. Lemmens, N. Ni, S. L. 

Bud'ko, and P. C. Canfield, Phys. Rev. B 78, 212503 

(2008). 

[37] L. Chauviere, Y. Gallais, M. Cazayous, A. Sacuto, and 
M.A. Measson, D. Colson, A. Forget, Phys. Rev. B 80, 
094504 (2009). 

[38] A. P. Litvinchuk, V. G. Hadjiev, M. N. Iliev, Bing Lv, A. 
M. Guloy, and C. W. Chu, Phys. Rev. B 78, 060503(R) 
(2008). 

[39] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 

[40] P. Giannozzi et al., http:/ /www. quantum-espresso. org" 
[41] T. Fukuda, A. Q. R. Baron, S. Shamoto, M. Ishikado, H. 

Nakamura, M. Machida, H. Uchiyama, S. Tsutsui, A. lyo, 

H. Kito, J. Mizuki, M. Aral, H. Eisaki, and H. Hosono, 

J. Phys. Soc. Jpn. 77, 103715 (2008). 
[42] T. Egami, B. Fine, D. J. Singh, D. Parshah, C. de la 

Cruz and P. Dai, Physica C, in press. 
[43] O. Delaire, M.S. Lucas, J. A. Muhoz, M. Kresch and B. 

Fultz, Phys. Rev. Letters 101, 105504 (2008) 
[44] O. Delaire, M. Kresch, J. A. Mufioz, M.S. Lucas, 

J.Y.Y. Lin, and B. Fultz, Phys. Rev. B 77, 214112 (2008). 



